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Abstract 

We consider a generalisation of the classical Lehmer problem about 
the parity distribution of an integer and its modular inverse. We use 
some known estimates of exponential sums to study a more general 
question of simultaneous distribution of the residues of any fixed num- 
ber of negative and positive powers of integers in prescribed arithmetic 
progressions. In particular, we improve and generalise a recent result 
of Y. Yi and W. Zhang. 

1 Introduction 

Given modulus q > 2, we denote by U q the set 

U q = {n : 1 < n < q, gcd(n, q) = 1}. 

In particular, j^lA q = (p(q), the Euler function. 

For n e U q we use n to denote the modular inverse of n, that is, the 
unique integer n &U q with nn = 1 (mod q) . 

The classical question of D. H. Lehmer (see [3l Problem F12]) is about the 
joint distribution of the parity of n and its modular inverse n e U q , defined 
by nn = 1 (mod q). 



1 



W. Zhang [T?l |T5] has shown that the Weil bound of Kloosterman sums, 
see Corollary 11.12], combined with some standard arguments, implies 
that if q is odd then n and its modular inverse n are of the same parity 
0.5<p(q) + O (y/2+o(i)) times for n e Uq 

This result has been extended in generalised in various directions, includ- 
ing its multidimensional analogues, see PQEl El El El [M HH H21 H31 [HI HS1 US] 
and references therein. 

In particular, it has been shown by Y. Yi and W. Zhang [IT] , that for 
any fixed integer k ^ 0, the smallest positive residue modulo q of n k and its 
modular inverse n k are of the same parity 0.5tp(q) + O (q 3 / i+ °^^ times for 
n E U q . 

Here we show that using the bound from [S] of exponential sums with 
sparse rational functions, one can get the same (as in the case k — 1) error 
term O (g 1 / 2 +°( 1 ) j for any fixed k and in fact obtain an asymptotic formula in 
a much more general case. Namely given an an integer s > 2 s-dimensional 
integer vectors 

k = (fa, . . . , k B ), ra = (mi,...,m s ), a = (a x , . . . , a s ), 

where mi, . . . , m s > 1, we denote by N q (m, a; k) the number of n e U q such 
that the smallest nonnegative residue of n k3 modulo q is congruent to aj 
modulo rrij for every j — 1, . . . , s. 

In particular, the result of Y. Yi and W. Zhang [11] can be reformulated 
as the asymptotic formula, 

N q ((2, 2), (0, 0); (k, k)) + N q ((2, 2), (1, 1); (fa k)) 

= l -y(q) + O (q^°V) , (1) 

which holds for any odd q. 

Here we give the following generalisation and improvement of (CE]). 

Theorem 1. For any fixed integer s > 2 and a vector k G Z s without zero 
components, uniformly over all vectors a, m e Z s with mi, . . . ,m s > 1 and 
an integer q > 1 with 

gcd(mi, q) = ... = gcd(m s , q) = 1, 

we have 

N q (m, a; k) = <p{q) + O (g^V'+od)) . 
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In particular, for s = 2, rri\ = m 2 = 2 and k\ = & 2 = Theorem [T] 
implies that the error term in the asymptotic formula ([TJ is O [q l l 2+o( - 1 ^ for 
every odd q > 1. 

Throughout the paper, the implied constants in the symbols 'O', and £ <C' 
may depend on the vector k. We recall that the notations U = 0(V) and 
V <C U are both equivalent to the assertion that the inequality \U\ < 
holds for some constant c > 0. 



2 Exponential Sums 



For an integer t we denote 

and recall that for hgZ, 
1 



e^(z) = exp(27U2:/£) 



e/(//u) 

-(e-l) /2<n<£/2 



1 ifu = (modf), 
otherwise, 



(2) 



(which follows immediately from the formula for the sum of a geometric 
progression) . 

We also recall that for any integers U > 1 and \x with < < £/2 we 
have 



u 



u=0 



<C min <( U, - — - 

m 



see [51 Bound (8.6)]. In particular 



E 

-(£-l)/2<IM<£/2 



u 



u=0 



< £logj 



(3) 



and 



E 

-(t-l)/2<jj,<e/2 



u=0 



< C/ + £log^ 



(4) 



Finally, as we have mentioned, our main tool is the following slight gen- 
eralisation of [H Theorem 1] which gives an estimate of exponential sums 
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with sparse rational functions. We note that in [HI Theorem 1] only the case 
of d = 1 has been considered, but the extension to the case of arbitrary d is 
st r aight f or war d . 

Lemma 2. For any fixed integer s > 2 and a vector k 6 I s without zero 
components and an integer q > 1, the bound 

neU q \l<j<s / 

holds, uniformly over all integers Ai, . . . , X s with 

gcd(Ai, . . . , X s ) = d. 

Proof. Let 



V 



Q = Y[Pi 

i=l 

be the prime number factorization of q and qi = q/p"% i = 1, . . . , v. We now 
define tj as the modular inverse of modulo pf\ that is, 

Uqt = (mod pf') and < U < pf% 

for i — 1, . . . , v. Using the multiplicative property of exponential sums with 
rational functions, see [5], Equation (12.21)] or [HI Lemma 6], we obtain 

H-iH v) =f[ E v f 4 - E ■ («) 

neU q \l<j<s J 1=1 nGW ttj \ 1<j'<s 

Furthermore, by [8], Lemma 5] we have 



L)) 



(6) 



for any prime power p Q with p a — >• oo and integers /ii, . . . , fi s with 

gcd(// x , . ..,fis,p) = 1. 
Combining ([5]) with (|SJ) and using that 

„ logg 

V < - — 

log log g 

(which follows from the obvious inequality v\ < q and the Stirling formula) 
we obtain the result. □ 
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2.1 Proof of Theorem [T] 



Without loss of generality we may assume that a has a nonnegative compo- 
nents satisfying < a, < mj, j = 1, . . . , s. 

Let us define Uj as the largest integers U with rrijU + aj < q, j = 1, . . . , s. 

Then N q (m, a; k) is equal to the number of solutions to the following 
system of congruences 

n k ' = rrijUj + ctj (mod q), n G U q , < Uj < Uj, j — 1, . . . , s. (7) 

Since gcd(m! . . . m s , q) = 1, for every j = 1, . . . , s we consider the mod- 
uluar inverse = mj of rrij modulo q, and also define bj e W g by the 
congruence 6j = djVj (mod g). Therefore, the system ([7]) is equivalent to the 
the following system of congruences 

rjn kj = + bj (mod g), n G W g , < Uj <Uj, j — 1, . . . , s. (8) 

Using (J2j) we write 



Changing the order of summation and then separating the main term 



iV 3 (m,a;k) = £ E - E 



n6W q 0<«i<!7i 0<u s <C/ s 




corresponding to Ai 



#U q Ui---Us _ y(g)Ui...U s 
q s q s 

. . = X s = 0, we obtain 



iVy(m,a; k) 
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where S* means that the term corresponding to Ai 
from the summation. Therefore, 



X s = is excluded 



N q (m, a; k) 



^{q)U x ...U s 



1 

< — 



E * 

-(9-l)/2<Ai,...,A s <9/2 



EM E J . 



new. 



n 

l<j<s 



E e *( A 



Now, for every divisor c? | q we collect together the terms with the same 
value of gcd(Ai, . . . , X s ) = d and then apply LemmaO obtaining the estimate 



N q (m, a; k) 



<p(q)U x ...U a 



< g -.+i-i/.+o(i)^ d i/. Efl 



q<q 



(9) 



where 



e n 



-(q-l)/2<\ 1 ,...,X s <q/2 l<j<s 
gcd(X 1 ,...,X s )=d 



E ma 



Writing \j = dfij, j = 1, . . . , s, and q = dqd, we derive 



e n 



-(?d-l)/2<Mi,-^ s <9ci/2 1<7'<s 
gcd(/xi,...,/i 3 )=l 



E e « 



Furthermore, we have 



where 



£<i < E ° d j 

i<j<s 



(10) 



e n 



-(gd-l)/2</ii,...,A» s <9d/2 l<j<s 0< Uj <Uj 



E e ^ (m] 



j = l,...,s. 
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We have, 



E 



'(9d-l)/2<M i <g d /2 



n E 



l<h<s -(q d -l)/2<ti h <q d /2 



0<u h <U h 



Applying ([3]) to the sum over fij and (TjJ to the other sums (and using the 
trivial estimate Uh < q), we obtain 

*dj < (q d logq d ) J] (C/ h + g d logg g ) < g^Qogg,) 5 = (TY +0 « 

l<h<s 



for every j = 1, . . . , s. Substituting this in f[10|) . and then recalling (J9J), we 
obtain 



iV 9 (m, a; k) - 



< g l-l/ S +0(l) ^ ^-1+1/* < gl-l/«+0(l) ^ 1 



By the well-known estimate on the divisor function, 



E 1 = 



see jH Theorem 317], we obtain 



N q (m, a; k) 



y{q)U x ...U s 



<q 



l-l/s+o(l) 



It remains to notice that Uj = q/rrij + 0(1), j = 1, . . . , s. Therefore 



Ut...U s 



q 



+ 0(q s ' 1 ] 



which concludes the proof. 
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